复杂动态环境下无人飞行器动态避障近似最优轨迹规划

针对复杂动态环境下无人飞行器的动态障碍规避问题，基于合理假设建立了无人飞行器和动态障碍的运动学模型，并综合考虑无人飞行器飞行过程中的终端约束、控制输入约束、安全避障约束等，以能量最少为性能指标构建动态避障问题数学描述。之后，针对终端约束和控制输入约束，依据优化模型预测静态规划算法(OMPSP)生成初始轨迹；针对动态避障问题的不等式约束，引入松弛变量并结合滑模变结构控制方法设计松弛变量动力学，实现对一个、多个或同时多个动态障碍的安全规避；最后，依据有限时间微分动态规划(RHDDP)算法进行轨迹优化，获得满足上述各种约束并能规避动态障碍的近似最优轨迹。     

A progressive approach to predict shot peening process parameters for forming integral panel of Al7050-T7451

In this paper, a progressive approach to predict the multiple shot peening process parameters for complex integral panel is proposed. Firstly, the invariable parameters in the forming process including shot size, mass flow, peening distance and peening angle are determined according to the empirical and machine type. Then, the optimal value of air pressure for the whole shot peening is selected by the experimental data. Finally, the feeding speed for every shot peening path is predicted by regression equation. The integral panel part with thickness from 2 mm to 5 mm and curvature radius from 3200 mm to 16000 mm is taken as a research object, and four experiments are conducted. In order to design specimens for acquiring the forming data, one experiment is conducted to compare the curvature radius of the plate and stringer-structural specimens, which were peened along the middle of the two stringers. The most striking finding of this experiment is that the outer shape error range is below 3.9%, so the plate specimens can be used in predicting feeding speed of the integral panel. The second experiment is performed and results show that when the coverage reaches the limit of 80%, the minimum feeding speed is 50 mm/s. By this feeding speed, the forming curvature radius of the specimens with different thickness from the third experiment is measured and compared with the research object, and the optimal air pressure is 0.15 MPa. Then, the plate specimens with thickness from 2 mm to 5 mm are peened in the fourth experiment, and the measured curvature radius data are used to calculate the feeding speed of different shot peening path by regressive analysis method. The algorithm is validated by forming a test part and the average deviation is 0.496 mm. It is shown that the approach can realize the forming of the integral panel precisely.     

基于双滤波器的高精度目标跟踪估计方法

针对连续推力的合作航天器，采用双重无迹卡尔曼滤波(DUKF)算法估计其状态和加速度。通过状态滤波器和参数滤波器的配合，提升滤波精度，完成运动状态和参数的估计，从而实现合作目标的运动轨迹跟踪。与合作航天器相比，非合作航天器存在大小未知、发生时刻未知的机动，无法获得加速度，且信息获取和运动状态的估计难度大。针对非合作航天器，基于简化的相对运动方程，结合天基平台获得目标的观测信息，采用两个扩展卡尔曼滤波(EKF)及基于半正焦弦的机动检测策略实现多未知脉冲机动的运动状态的估计。仿真结果表明：相比于无迹卡尔曼滤波（UKF），DUKF在对合作航天器的状态和加速度估计方面具有更快的收敛速度和更高的滤波精度；对于存在未知机动的非合作航天器，通过对比验证机动检测策略与滤波器切换策略相结合的方法的有效性，该方法能够检测到多次机动并且减少误判。     

航天嵌入式软件数组越界缺陷特征研究

根据统计,数组越界是航天嵌入式软件开发过程中出现最多且最容易被遗漏的缺陷类型之一.目前自动化检测数组越界多基于抽象解释、符号执行、程序模型检验等方法,这些方法在误报、漏报、可扩展性等方面的表现依赖于软件及缺陷特征.分析了近三年航天嵌入式软件第三方测试中发现的94个数组越界问题,从缺陷模式和缺陷表现形式两方面分析得出10项航天嵌入式软件数组越界缺陷特征,并提出对设计具体检测方法关键的若干启示.进一步基于这些特征和启示探讨了数组越界检测算法针对中断驱动型程序的改进方向.     

星载激光通信终端在轨机动对温度影响

对激光通信终端在轨瞬态温度变化开展了仿真，以期研究在轨机动的影响、热分析和热试验时机动模式的简化模拟。通过合理地分析与简化，建立终端的轨道热分析模型，准确模拟在轨机动，根据典型机动模式、外热流和涂层退化等因素设计了计算工况，得出了终端在轨运行过程中的温度场随时间和姿态变化规律。结果表明：不同机动模式下的温度变化存在差异，最大可达23.0℃，采用固定姿态或一维转动的简化热分析或热试验不能准确模拟实际飞行温度，甚至不能部分替代二维转动热分析或热试验；光学天线和反射镜的温度控制是制约终端工作的瓶颈，为终端设计合适的避光机动策略，可大幅度提高温度稳定度和均匀度。研究结果可以为光机电设备在轨机动策略的设计和改进提供参考。     

Topology optimization in lightweight design of a 3D-printed flapping-wing micro aerial vehicle

Topology optimization is an effective method to obtain a lightweight structure that meets the requirements of structural strength. Whether the optimization results meet the actual needs mainly depends on the accuracy of the material properties and the boundary conditions, especially for a tiny Flapping-wing Micro Aerial Vehicle (FMAV) transmission system manufactured by 3D printing. In this paper, experimental and numerical computation efforts were undertaken to gain a reliable topology optimization method for the bottom of the transmission system. First, the constitutive behavior of the ultraviolet (UV) curable resin used in fabrication was evaluated. Second, a numerical computation model describing further verified via experiments. Topology optimization modeling considering nonlinear factors, e.g. contact, friction and collision, was presented, and the optimization results were verified by both dynamic simulation and experiments. Finally, detailed discussions on different load cases and constraints were presented to clarify their effect on the optimization. Our methods and results presented in this paper may shed light on the lightweight design of a FMAV.     

轨控大干扰下的系统角动量管理

摘要： 针对利用角动量交换装置系统实现轨控期间姿态维持的情况，研究当轨控推力器存在大力矩扰动时系统角动量管理方法.结合对地运行航天器轨道运动特性，分析轨控干扰力矩的积累角动量变化规律；考虑角动量管理装置系统的角动量存贮容量，基于角动量积累规律提出了一种结合偏置角动量建立与对称分布轨控位置选择的系统角动量管理方法，实现了对称轨控位置积累角动量互相抵消的自平衡效果，极大程度地提高轨控效率且避免系统角动量饱和现象.所提出的方法适用于喷气欠驱动控制系统，方法的有效性通过了数学仿真验证及在轨型号应用.     

Accurate approximation of in-ecliptic trajectories for E-sail with constant pitch angle

Propellantless continuous-thrust propulsion systems, such as electric solar wind sails, may be successfully used for new space missions, especially those requiring high-energy orbit transfers. When the mass-to-thrust ratio is sufficiently large, the spacecraft trajectory is characterized by long flight times with a number of revolutions around the Sun. The corresponding mission analysis, especially when addressed within an optimal context, requires a significant amount of simulation effort. Analytical trajectories are therefore useful aids in a preliminary phase of mission design, even though exact solution are very difficult to obtain. The aim of this paper is to present an accurate, analytical, approximation of the spacecraft trajectory generated by an electric solar wind sail with a constant pitch angle, using the latest mathematical model of the thrust vector. Assuming a heliocentric circular parking orbit and a two-dimensional scenario, the simulation results show that the proposed equations are able to accurately describe the actual spacecraft trajectory for a long time interval when the propulsive acceleration magnitude is sufficiently small.     

拦截高速目标的比例与反比例导引捕获区分析

针对拦截高速目标的作战特点，分析了比例导引(PN)与反比例导引(RPN)的捕获区。首先，通过分析拦截弹与目标的相对运动关系，推导得到了顺轨和逆轨的零控拦截条件，此条件由目标和拦截弹的速度前置角以及二者速度比确定；其次，以拦截弹和目标速度前置角为坐标系，推导得到了PN以及RPN捕获区以及各自导航比设置范围。PN的捕获区由逆轨零控拦截条件以及与其相切且斜率为1/(N-1)的两条直线构成，RPN的捕获区由顺轨零控拦截条件以及与其相切且斜率为1/(-N-1)的两条直线构成；然后，利用函数对称性将PN与RPN捕获区转换到同一坐标区间，得到了相同条件下RPN捕获区要大于PN捕获区的结论；最后，开展了四种情形下的仿真，验证了本文捕获区分析的合理性及有效性。     

不同密度弹丸对水冰的超高速撞击成坑实验

为了研究冰冻天体表面撞击坑的形成与演化，开展了水冰的超高速撞击成坑实验。使用二级轻气炮发射1.0 mm直径的球形弹丸，以3 km/s、5 km/s和7 km/s速度对圆柱状冰块进行撞击。弹丸材料包括聚碳酸酯和不锈钢两种，冰块温度为253 K。实验观察到了不同弹丸和不同速度条件下，冰块中撞击坑的形貌特征。对撞击坑直径、深度和剖面形状进行了测量，并与文献中铝弹丸对水冰的撞击坑进行了比较分析。获得了水冰撞击坑特征随撞击参数的变化规律，结果表明：撞击坑直径和深度的主导机制不同，坑深主要由弹丸侵彻作用形成，而坑径主要由冰块的剥落所致；坑深比坑径具有更强的对于弹丸密度的依赖性，高密度弹丸撞击坑直径具有比低密度弹丸更强的对于撞击速度的依赖性；撞击坑体积与撞击能量成正比，高密度弹丸形成的撞击坑直径表现出“能量缩比”行为，而低密度弹丸形成的撞击坑直径表现出“动量缩比”行为。